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PREFACE 


This  report  is  generated  under  Contract  DAAK70-79-C-0158  for  the  U.S.  Army 
Engineer  Topographic  Laboratories,  Fort  Belvoir,  Virginia  22060,  by 
Autometric,  Inc.,  Falls  Church,  Virginia,  and  submitted  as  ETL-0306.  The 
Contract  Officer's  Representatives  were  Mr.  John  Armlstead  and  Mr.  Eugene 


1.0 


INTRODUCTION 


The  objective  of  this  effort  is  to  mathematically  model  the  two 
media  refraction  process  and  incorporate  this  model  into  the  existing  MUSAT 
aerotrlangulatlon  programs.  This  will  permit  these  revised  programs  to 
accomodate  photographic  image  points  that  appear  in  shallow  water  as  well  as 
those  which  appear  above  the  surface. 

Classic  hydrographic  mapping  in  shoal  waters  using  surface 
vessels  is  hazardous,  costly,  and  very  time  consuming.  Precision  three- 
dimensional  photogrammetric  techniques  offer  a  basic  alternative  which 
minimizes  these  drawbacks*  Single  stereophotogrammetrlc  models  are  currently 
being  oriented  in  analog  plotters  for  the  compilation  of  depth  curves  and 
bathymetric  features.  However,  this  can  only  be  done  where  sufficient  control 
exists  to  scale  and  orient  the  model.  Photogrammetric  aerotrlangulatlon 
techniques  currently  in  use  would  appear  to  satisfy  the  need  for  control 
denslflcation  in  bathymetric  models.  With  the  mathematical  model  in  place, 
large  blocks  of  overlapping  aerial  photographs  can  be  used  to  map  extensive 
offshore  shoal  areas. 

In  considering  the  development  of  the  math  model  to  correct  for 
refraction,  the  following  assumptions  have  been  made: 

a.  The  elevation  of  the  water  surface  is  known  from  tidal 

data,  for  each  exposure. 

b.  The  wave  induced  refraction  effects  will  result  in  plate 

residuals  and  ground  point  errors  which  are  uncorrectable 
in  this  effort. 

c.  Light  rays  reflected  from  a  bottom  feature  travel  in  a 

straight  line  while  in  the  water. 

d.  As  a  result  of  these  assumptions,  there  will  be  some 


discrepancies  which  cannot  be  accounted  for.  Within 
reason,  these  will  be  included  in  the  normal  plate 
residuals. 

This  final  report  presents  the  results  of  the  modifications  to 
the  MUSAT  aerotrlangulatlon  programs  to  accommodate  bathymetric  image  points, 
as  contracted  by  USAMERADCOM.  Two  versions  of  the  MUSAT  programs  were 
modified,  MUSAT-IIA,  and  MUSAT-IIIA.  These  two  programs  use  different  theory 
in  their  solution  algorithms  and  therefore  warrant  different  approaches  to  the 
multi-media  refraction  problem  incurred  by  bathymetric  points. 

The  MUSAT-IIA  software  package  utilizes  the  theory  of  modified 
geodetic  restraint.  Under  this  theory,  only  camera  station  parameters 
(attitude  and  position)  are  considered  unknown.  The  collinearlty  condition 
equations  are  then  written  for  all  control  points  whether  fully  known  or 
not.  Coplanarity  and  scale  restraint  are  then  applied  to  all  pass  points. 
Equal  elevation  and  known  airbase  conditions  can  also  be  enforced. 

The  MUSAT-IIIA  program  performs  a  simultaneous  block  solution 
using  all  points.  Collinearlty  condition  equations  are  written  for  all  points, 
whether  fully  known  or  not.  The  system  of  equations  is  solved  simultaneously 
for  all  point  coordinates.  Consequently,  the  water  refraction  corrections  must 
be  applied  to  the  observed  image  coordinates  for  every  iteration  of  the 
simultaneous  solution. 

The  remainder  of  the  report  discusses  the  changes  to  the  input 
requirements  and  output  formats  for  the  bathymetric  solution  for  both  MUSAT- 
IIA  and  MUSAT-IIIA.  The  modifications  to  the  existing  AS-11AM  software  to 
allow  compilation  of  subsurface  Information  are  also  presented. 


All  programming  changes  and  additions  to  the  MUSAT  programs  are 
written  in  FORTRAN  for  Implementation  on  DMA/HTC  Uni vac  1108  under  the  EXEC  8 
operating  system.  Program  changes  to  the  AS-11AM  software  are  written  in 
FORTRAN  for  implementation  on  the  M0DC0MP  11/25  under  the  MAX  III  operating 

system. 


2.0 


Theory  of  Multi-Media  Refraction 


I 


From  a  photogrammetric  standpoint,  the  effect  of  refraction  due 
to  the  water  column  above  a  submerged  point,  is  a  displacement  of  the  point's 
image  on  a  photograph.  The  submerged  point's  image  is  displaced  radially 
outward  from  the  principal  point.  The  magnitude  of  the  displacement  is  a 
function  of  the  spatial  relationship  existing  between:  the  exposure  station 

attitude  and  position,  the  submerged  point,  and  the  height  of  the  wat 
surface,  at  the  time  of  exposure. 


The  law  may  be  stated  as: 


sin  6 


sin  I 
n 


where : 


n 

I 

6 


Refractive  index  at  the  interface 
Angle  of  incidence  of  the  image  ray 
Angle  of  refraction  of  the  image  ray 


It  Is  appropriate  at  this  point  to  stress  that  no  attempt  was 
made  to  model  the  effects  of  wave  action  and  that  all  refractive  corrections 
assume  a  smooth  surface.  Figure  2.2  illustrates  the  effect  of  water 
refraction  in  a  stereo  model. 


tide  height  photo  1 


Focal  lengths  for  each  photograph 


At  this  point,  we  can  begin  the  derivation  of  the  mathematical 
models  for  refraction.  As  noted  earlier,  due  to  the  differences  between 
MUSAT-IIA  and  MUSAT-IIIA,  different  models  are  required  for  each  version. 
Each  will  be  derived  independently  to  retain  clarity. 


3.0 


MUSAT-IIA 


3.1  Derivation  of  Mathematical  Model  /or  Refraction  Correction 


The  first  step  in  the  algorithm  concerns  the  establishment  of  a 
local  vertical  coordinate  system  sufficiently  near  the  submerged  points' 


position  so  as  to  not  be  affected  by 
convention  is  established: 

Xi,  Y±,  Z± 

Mi 

fi 

x'i*  y'i 

xp«»  YP’»  ZP’ 


earth  curvature.  The  following  notation 

-  Camera  station  "i"  geocentric 

position 

Camera  to  geocentric  rotation  matrix 
for  exposure  station  ”i" 

-  Focal  length  of  camera  "i" 

-  Image  coordinates  of  an  apparent 
submerged  point  on  photograph  ”i” 

-  Geocentric  position  of  point  P' 


where : 


The  collinearity  condition  is: 


(3.1) 


Scale  factor 


-  Object  space  image  ray  corresponding 
to  the  apparent  position  of  the 
submerged  point. 


(It  is  assumed  that  the  focal  length  carries  the  proper  algebraic  sign) 

By  intersecting  the  Image  rays  from  the  first  two  exposure 


stations,  we  can  define  an  appropriate  origin  for  a  local  vertical  coordinate 
system.  The  mathematics  for  this  intersection  are: 


(3.2) 


(D  x  B)  RI 
k.  -  - - 

1  (B  .  B) 

(D  x  B) 

k  -  - L 

2  (B  .  B) 


(3.3) 


RJ  x  R^ 


X2-Xl 
Y2  -  Y1 
Z2-Zl 


+  k  h  +  Y„ 

11  4 


+  k2*2 


(3.4) 


(3.5) 


(3.6) 


We  can  now  convert  the  geocentric  coordinates  of  the  point  P' 
into  geographic  coordinates  (4>pI,  Xpl,  hpt).  Also,  the  local-to-geocentric 
orientation  matrix  (  M  )  can  be  formed.  This  matrix  takes  the  form: 

-sin  X  -sin  *  cos  X  cos  6  cos  X 

p«  p«  p»  pi  p 

M  -  cos  X  -sin  <j>  sin  X  cos  4  sin  X  (3.7) 


0  cos  d>  sin  d> 

L  P’  P’  J 

This  orientation  matrix  defines  a  local  system  in  which  the  "Z" 
axis  is  the  local  vertical.  Making  P'  the  origin  of  the  local  system,  the 
transformation  from  local  to  geocentric  is  now: 


where : 


Y- V 

z  -  zpl 


x,  Y,  Z 

*.  y,  * 


x 

M  y 


(3.8) 
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Geocentric  coordinates  of  a  point 
Local  coordinates  of  the  point 


The  next  task  Is  to  compute  the  point  of  Intersection  of  all 
Image  rays  for  the  submerged  point  (not  just  those  on  the  first  two 
photographs)  with  the  water  surface.  Using  equations  (3.8)  and  (3.1)  the 
expression  for  the  conversion  of  an  image  ray  to  the  local  system  is: 


Yi  +  kiMi  y'i 


Yp,  +  M 


(3.9) 


V  V  T  xi 

Yi"  yp»  ■  ki  * \  yl 
zi-  V  fi 


(3.10) 


Further  substitution  yields  a  compact  solution  for  the  object 
space  image  ray  at  exposure  station  "i"  in  the  local  coordinate  system: 


y  -  y±  +  kjfi 


yi  +  ki  ri 


(3.11) 


where: 


xi’  yi’  Zi 


-  Local  coordinates  of  exposure 
station  ”i" 


-  Direction  numbers  of  object  space 
image  ray  "i”  in  the  local  system 


x,  y»  * 


-  Terminus  of  the  object  space  image 
ray  in  the  local  system 


The  intersection  of  object  space  ray  ”i”  with  the  water  surface 


can  now  be  defined  by  assigning  "z"  a  value  equal  to  the  tide  elevation 


f 

i 


associated  with  photograph  "i",  in  the  local  coordinate  system.  Therefore, 


Geographic  height  of  the  tidal  surface. 


Equations  (3.11)  and  (3.12)  are  combined  to  yield  a  solution  for  the  water 
surface  Intersection: 


zt)  (r1/r3) 


*1  + 


z±)  (r2/r3) 


(3.13) 


where: 


Elements  of  the  direction  number 
♦ 

vector  r^ 

Local  coordinates  of  the  water 
surface  Intersection 


The  Law  of  Refraction  is  now  applied  to  each  object  Bpace  ray  by 
first  computing  the  angle  of  incidence  of  the  ray  and  then  determining  the 
angle  of  refraction.  This  Is  done  using  the  following  expressions: 


tan(  I  )  -  — - 

Zi  "  (hti"  bp'* 

5i  ■  «;i  -  ;.l)2  +  ‘i-i  -  \>2>/z 


sin 


-1 


sin(I^) 


(3.14) 


(3.15) 


(3.16) 


The  object  space  ray  under  the  water  surface  can  now  be 


determined  for  each  Incident  ray  by: 


(h  -  h  ,)  tan(0  )  (x  -  x  ) 
i  l. 


(ht1  “  hp,)  tan(0i>  <y„  "  yi) 


(3.17) 


where: 


-  v> 


-  Underwater  object  space  ray  in  the 
local  coordinate  system 


We  now  have  the  situation  depicted  in  Figure  3.1.  With  "m"  image 
points  we  can  form  "m”  underwater  object  space  rays  **  u^  ”.  A  new  estimate  of 
the  submerged  point's  position,  P",  can  be  determined  by  the  intersection  of 
all  the  rays.  The  equations  required  take  the  form  of  the  collinearity 
condition: 


x_„  -  x 
P  w. 


ypM  ”  V 


ki  *  “i 


(3.18) 


z_„  -  z 
P  w. 


A  quick,  and  sufficiently  rigorous  intersection  can  be  performed  by 


rearranging  equation  (3.18)  to  become: 


where: 


m  (ui^u3\ 

yP"  ”  (u2/u3)i  <V> 


Wi  <2»l> 

Wi  (\> 


(3.19) 


Uj  f  U2  9  Uj 


-  Elements  of  the  underwater  object 
•pace  ray  ”  u.  ” 
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The  result  Is  the  following  system  of  equations  for  all  object  space  rays  for 
the  submerged  point: 


1  0  -(Uj/u^j 

o  i  -oyupj 


\  '  <ul/u3)l<*w1) 
yw  ~  ^u2^u3^1^zw  ^ 


(3.20) 


1  0  -(u1/u3)m 

0  1  -Cu2/u3)m 


V  -  (ul/u3)m(zw  > 

1  m 

v  -  'V-sV*,  > 

1  m 


The  solution  of  this  system  of  equations  yields  a  new  estimate  of  the  position 
of  the  submerged  point  P",  in  the  local  system.  The  new  coordinates  may  be 
converted  to  geocentric  by: 


Xp..  Xp, 


V  -  YP.  +  m  yp. 


(3.21) 


The  final  step  of  the  algorithm  is  to  determine  a  new  set  of 
Image  coordinates  for  each  image  of  the  submerged  points.  The  Image 
coordinates  are  computed  by  rearranging  the  colllnearlty  condition  equations: 


(3.22) 


(3.23) 
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-  Bn(Xp..  -  \)  +  ®12(V  -  Yi>  +  ®13<ZP"  "  zi )  (3.24) 

Oj  ■  ®2i^p"  ”  ^  ®22^P”  ”  T^)  ^  m23^ZP”  ”  Zi ^  (3.25) 

*  ®3^(Xp'>  _  3C^)  +  ®32^YP"  ”  Y^ )  +  ®33(zp"  "  Zi^  (3.26) 

where: 

nijj  -  m^  -  Elements  of  the  camera  to  geocentric 

rotation  matrix 

x^,  y^  -  New  image  coordinates  of  the 

submerged  point  based  on  the  updated 
(true)  position  of  the  submerged 
point 


It  can  be  shown  that  as  the  difference  in  position  of  the  points 
P'  and  P"  becomes  neglibly  small,  that  P”  actually  approaches  the  true 
position  of  the  submerged  point,  P.  This  criteria  will  be  used  to  test  for 
convergence  of  the  water  point  correction  solution. 
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3.2 


Description  of  New  Subroutine 


A  new  subroutine  has  been  incorporated  into  the  MUSAT-IIA 
aerotriangulation  software  package.  It  corrects  the  observed  image 
coordinates  of  an  apparent  submerged  point  to  reflect  its  true  position. 
Subroutine  WATPT  is  engaged  only  when  a  submerged  point  is  encountered  in  the 
input  data  for  a  triangulation  run.  Should  there  be  no  submerged  points,  its 
inclusion  will  in  no  way  affect  the  solution  of  a  standard  job. 


3.2.1 


SUBROUTINE  WATPT 


SUBROUTINE: 

DATE  WRITTEN : 
LANGUAGE/MACHINE : 
DESCRIPTION: 

CALL  LINE: 

INPUT: 

OUTPUT: 

METHOD: 


WATPT 

November  1980 

FORTRAN  /  UNIVAC  1108  Exec  8 

This  routine  corrects  the  observed  image 
coordinates  of  a  bathymetric  point  to  reflect  it's  true 
position. 

Call  WATPT 

Thru  ITAPE3 : 

GNDINF  -  Image  point  records 


Thru  ITAPE2 : 

GNDINF  -  Image  point  records  with  new 

estimates  for  the  image  cordinates 

As  stated  earlier,  the  image  coordinates 
corrections  are  made  to  a  pre-triangulated  data  set,  and 
therefore  estimates  exist  for  all  camera  and  ground  point 
parameters.  The  routine  cycles  through  each  physical  point 
record  and  skips  it  if  no  submerged  point  is  found.  Upon 
locating  a  point  record  for  a  submerged  point  the  following 
operations  take  place. 


First,  a  local  vertical  coordinate  system  is 
established  at  the  intersection  of  the  image  rays  of  the 
submerged  point  from  the  first  two  photographs  on  which  the 
point  appears.  Then,  the  image  rays  from  all  photographs 
on  which  the  point  appears  are  intersected  to  produce  a 
better  estimate  of  the  point's  position.  All  relevant 
parameters  are  now  transformed  into  the  local  vertical 
system,  including  the  tidal  elevation.  The  resulting 
coordinates  for  the  submerged  point  are  transformed  into 
the  geocentric  coordinate  system. 


At  this  point,  updated  image  coordinates  are 
computed  for  each  image  of  the  submerged  point,  reflecting 
the  effect  of  the  change  in  the  point's  position  on  its 
corresponding  images.  These  corrected  image  coordinates 
and  submerged  point  coordinates  are  written  out  to  the 
physical  point  record.  This  process  is  repeated  for  each 
submerged  point  encountered  as  the  routine  cycles  through 
the  physical  point  records. 
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LIMITATIONS:  NONE 

COMMON  BLOCKS:  /WATER/  /TAPES/  /CAMDAT/ 

/AXES/  /TOTAL/ 

OPTIONS:  NONE 

SUBROUTINES 

REFERENCED:  Function  BIGM,  MPYAB,  SUBMAT,  ADDMAT,  CROSS,  DOT,  XYZPLH, 

ROTFRM,  INVRT 
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3.3 


DICTIONARY  OF  NEW  COMMON  BLOCKS  -  MUSAT-IIA 


VARIABLE 


IWTPT 


REFIDX 


COMMON  /WATER/  IWTPT,  REFIDX 

TYPE  DESCRIPTION 

Integer  Bathymetric  points  flag 

0  *  No  bathymetric  points  this  job 
1  -  Bathymetric  points  in  this  job 

Real  Index  of  refraction  of  the  water  body 

Nominal  values: 

n  »  1.3333  Fresh  water 


n  “  1.3334  Sea  water 


4.0 


MU SAT- I I I A 


4.1  Derivation  of  Mathematical  Model  for  Refraction  Correction 

The  mathematical  algorithm  to  compensate  for  the  refracted  images 
of  underwater  points  must  address  two  problems.  The  first  defining  a  local 
level  surface  in  each  of  the  various  coordinate  systems  utilized  by 
MUSAT-IIIA.  The  second  problem  is  modifying  the  mathematics  of  the  solution 
for  water  penetration. 

The  MUSAT-IIIA  program  contains  preprocessors  for  blunder  editing 
and  for  obtaining  initial  approximation  for  the  solution.  It  has  been  assumed 
that  these  approximations  are  accurate  enough  to  perform  a  rigorous  solution 
taking  into  consideration  the  effect  of  water  penetration. 

The  first  mathematical  problem  addressed  Is  the  definition  of  the 
local  level  surface  for  each  underwater  point.  Since  the  position  of  each 
point  Is  not  known,  this  is  not  a  simple  matter.  In  the  preprocessing  stage, 
the  initial  approximation  for  the  position  is  determined.  Once  the 
preprocessing  is  complete,  a  coordinate  transformation  is  set  up  for  each 
underwater  point  and  stored  In  the  corresponding  point  file.  This 
transformation  is  a  conversion  between  the  triangulation  coordinate  system  and 
a  local  vertical  coordinate  system  at  each  point. 

Because  of  the  earth  curvature  problems ,  and  the  problem  of 
defining  the  local  vertical,  It  is  assumed  that  all  water  penetration 
triangulation  will  be  performed  in  an  earth  centered  coordinate  system,  and 
that  the  actual  reduction  will  be  In  some  cartesian  system,  either  geocentric 
or  local  tangent. 

As  the  points  are  processed,  the  positions  of  these  points  change 


with  each  iteration  of  the  solution.  For  a  completely  rigorous  solution, 
during  the  point  processing,  the  transformation  matrix  for  each  point  is 
recomputed  when  the  point  position  is  updated  for  each  iteration. 

In  order  to  take  advantage  of  the  geometry  of  light  refraction  in 
water,  computation  oust  be  made  in  a  local  vertical  coordinate  system 
established  for  each  point.  The  transformation  to  these  coordinate  systems, 
as  discussed  above,  will  be  used  to  convert  all  coordinates  to  a  local 
system.  In  the  geometry  of  this  local  system,  the  modified  collinearity 
equations  can  be  established.  The  partial  derivatives  of  the  collinearity 
equations  will  then  be  converted  to  the  triangulation  coordinate  system. 

The  following  figure  shows  the  geometry  of  the  refracted  light 
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It  is  assumed  for  the  solution  that  the  elevation  of  the  water,  Z  Is  a  known 

w 

quantity. 

The  colllnearity  condition  does  not  hold  for  a  submerged  point 
due  to  the  refraction  Imparted  by  the  water  column  above  the  point.  The 
general  colllnearity  equations  will  be  modified  to  reflect  the  geometry  in 
Figure  4.1.  The  colllnearity  equations  are  written  as: 


X  -  X- 
w  C 

Y  -  Y  _ 
w  C 

Z  -  Z. 
w  C 


-  k  R  y  -  y 


(4.1) 


where , 


V  yP 


-  Image  coordinates  of  a  point 

-  Principal  point  offset 

-  Focal  length 

-  Photo  to  ground  rotation  matrix 

-  Scale  factor 

-  Coordinates  of  image  ray/water 
surface  intersection 


XC  ,YC  ,ZC 


-  Camera  station  coordinates 


However,  Xy.  Y^^  are  functions  of  Xj,,  Yg,  Z^,  and  Xq,  Yq,  Zq.  By  differen¬ 
tiating  the  functions  of  Xw,Yw  with  respect  to  the  camera  station  and  ground 
point  coordinates;  the  variables  Xw»Yy  can  be  eliminated  from  the  equations. 
The  partial  derivatives  of  the  projective  equations  will  now  be  in  terms  of 


X„,  Y„,  Z„,  X-,  Y„,  Z„,  as  they  should  be. 


First,  we  define  an  axis  Q  such  that  the  light  ray  lies  in  the  QZ 
plane  with  the  origin  at  the  water  surface.  Viewing  the  light  ray  from  above 
the  XT  plane,  the  ray  would  appear  as  in  Figure  4.2.  From  Figure  4.1  and  4.2 
we  define  the  following: 


Figure  4.2  Top  View  of  a  Bathymetric  Image  Ray. 
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Qc  -  ((Xg  -  X„)2  +  <Yl  - 

°G  "  <(Xw  -  *g)2  +  (Yw  - 

D  -  «Xc  -  Xg)2  +  <Yl  - 

It  can  be  seen  that, 

D  -  Qc  +  0G 

and 

QG  “  D  -  Qc 


V2^2 

(4.2) 

v2//2 

(4.3) 

CM 

(4.4) 

- 

(4.5) 

(4.6) 

From  Figure  4.1, 


Sin  I 


v2 


+  (Z  - 
w 


v2> 


(4.7) 


Sin  6 


<V  +  (ZG  -  V  > 


D~QC 

( (D  —  Og)2  +  (ZG  -  Zw)2)  X/2 


(4.8) 


Snell's  Law  of  Refraction  states: 


Sin  I  ■  n  Sin  0 


(4.9) 


where , 

n 


Index  of  refraction 


The  function  for  combining  all  the  parameters  is  then  : 
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■In 


<0c2  +  (z. -V2)i/2 


-  n  sin 


(D  "  °C) 

((D  -  Qc)2  +  (ZG  -  Zw)2)1;2 


(4.10) 


The  function  F  can  now  be  evaluated  by  iteration  for  the  current  values  for 

V  V 

Let  us  define  the  colllnearity  equations  as  follows: 


(x  -  xp)  +  f 


R11  (Xw  -  V  +  R12  (Yw  -  V  *  R13  Uw  -  ZC) 
R31  (Xw  *  V  +  R32  (Yw  r  V  +  R33  (Zw  '  V 


0 

(4.11) 


F 

y 


(y  -  yp)  +  f 


R2i  Uw  -  xc)  +  R22  (Yw  -  V  +  R23  (2w  -  zc) 
R31  (XW  -  V  +  R32  (YW  -  V  +  R33  (ZW  "  ZC> 


0 

(4.12) 


or. 


(x  -  xp)  +  f  ^ 


(4.13) 


Fy  -  (y  -  yp)  +  f  w  “  0 


(4.14) 


In  order  to  remove  1^,  Yw  from  the  partlals,  the  chain  rule  oust  be  applied  in 
computing  the  following  partlals: 


-24- 


Differentiating  the  refraction  function,  F,  using  the  chain  rule,  we  have 


3F 

3F 

30c 

3F 

3D 

" 

3xc 

+  3D 

3F 

3F 

3Qc 

3F 

3D 

‘ 

3Yc 

+  3D 

3F 

3F 

30c 

3F 

i 

3D 

3zc 

3D 

Similarly,  we  have: 


3F 

3F 

3QC 

« 

3F 

3D 

3X~ 

w 

"  Wc 

3X 

w 

T 

75” 

75T 

w 

3F 

3F 

3°c 

i 

3F 

3D 

3Y~ 

w 

" 

3Y 

w 

■r 

TIT 

3Y” 

w 

3F 

3F 

3QC 

* 

3F 

3D 

Tz~ 

30l 

3Z. 

*r 

75” 

3Z 

(4 


(4 


(4 


We  can  combine  these  partlals  to  form  the  required  expressions.  Specifically, 


3F 

3F 

3F 

3X 

w 

3xc 

3xw 

3Yc 

3X 

w 

3zc 

3xc 

3F 

3Yc  " 

3F 

3Zc  " 

3F 

3X 

3X 

3X 

w 

w 

w 

3F 

3F 

3F 

3X 

3 

3  X 

3  Y 

3X 

3Z„ 

w 

G 

w 

G 

w 

G 

3xg  " 

3F 

3YG  ’ 

3F 

3Zr  " 

3  F 

3X 

w 

la 

3  X 
w 

la 

3X 

w 

Similar  construction  results  in  the  partials  of  Yw,  and  Zw  with  respect  to  the 
camera  and  ground  point  coordinates. 


Now,  taking  the  partial  derivative  of  the  collinearity  equations 
with  respect  to  camera  station  position,  for  example,  and  adding  the 
contribution  of  the  refraction  function,  the  full  partial  becomes: 


x 

W 


V 


(4.23) 


31 


) 


+ 


3F  3Y  3F  3Z 

__Z  _ S  +  _Z  -JL 

3  Y  3  X  3  Z  3  X_ 

w  C  w  C 


(4.24) 
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In  the  manner  shown  above,  the  discrepancy  vector  and  the  partial 
derivatives  will  be  modified  for  all  points  which  are  under  the  water  on  each 
photograph.  This  will  be  accomplished  by  a  modification  of  the  routines  which 
compute  these  discrepancies  and  partlals,  which  are  used  universally 
throughout  the  program.  The  solution  and  its  attendant  error  propagation  of 
the  errors  through  various  output  coordinate  systems  will  remain  valid. 
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4.2  Description  of  New  Subroutines 

Three  new  subroutines  have  been  incorporated  Into  the  MUSAT-IIIA 
software  package.  These  new  routines  deal  specifically  with  the  modifications 
made  to  the  existing  software  to  accommodate  bathymetric  image  points.  Should 
there  be  no  requirement  for  this  capability  in  any  specified  MUSAT-IIIA  run, 
their  Inclusion  will  in  no  way  affect  the  solution  of  the  standard 


aerotrlangulatlon  problem 


4.2.1 


SUBROUTINE  UNWTR 


SUBROUTINE: 

DATE  WRITTEN 
WRITTEN  BY 
LANGUAGE /MACHINE : 
DESCRIPTION: 

CALL  LINE: 

INPUT: 


OUTPUT: 


CALLED  BY: 
METHOD: 


UNWTR 
June  1980 

L.  Peters  Autometric.  Inc. 

FORTRAN  /  UNIVAC  1108  Exec  8 

This  routine  computes  the  coefficients  of  the 
colllnearity  condition  equations  for  the  submerged 
points.  In  this  sense,  the  routine  is  functionally  similar 
to  the  subroutine  PTCOND  in  the  existing  MUSAT-IIIA 
software. 

CALL  UNWTR  (  ICOL,  WATER,  XYZ,  ISW  ) 

Thru  the  argument  list  - 

ICOL  -  Column  pointer  to  an  image  in  the 

physical  record 

XYZ  -  Ground  point  coordinates  in  the 

triangulation  system 

ISW  -  Operational  switch  governing  the  use 

of  the  routine 

Thru  the  argument  list  - 

XWATER(3)  -  Image  ray /water  surface  intersection 

coordinates 

Thru  COMMON  - 


ATTCOE 

-  Attitude 
condition 

coefficients 

equations 

of 

the 

POSCOE 

-  Position 
condition 

coefficients 

equations 

of 

the 

GRNCOE 

-  Ground  coefficients  of  the 
equations 

condition 

CONCOE 

-  Constant 
condition 

coefficients 

equations 

of 

the 

LNK6C1,  LNK6C2 


The  first  task  addressed  is  to  check  the  input 
coordinate  system  and  perform  the  necessary  coordinate 
transformations  to  reduce  it  to  a  geocentric  coordinate 
system.  This  procedure  Includes  transforming  both  the 
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LIMITATIONS: 
COMMON  BLOCKS: 

OPTIONS: 

SUBROUTINES 


coordinates  and  the  orientation  matrices  as  necessary. 
These  coordinates  are  transformed  into  a  local  vertical 
system  with  its  origin  at  the  initial  approximation  for  the 
submerged  point's  position  as  generated  by  the  MUSAT 
preprocessors.  We  now  have  all  the  appropriate  parameters, 
l.e.,  camera  station  parameters,  ground  point  parameters, 
and  associated  orientation  matrices  in  this  local  vertical 
system.  At  this  point  the  value  for  the  geographic  height 
of  the  tidal  elevation  is  also  transformed  to  this  system. 


The  next  task  is  to  compute  the  coordinates  of 
the  object  space  image  ray  -  tidal  surface  intersection 
point.  This  is  performed  by  invoking  the  second  new 
subroutine,  WATPNT,  which  will  be  discussed  in  detail  in 
succeeding  sections.  The  results  of  this  routine  are  the 
X,Y  coordinates  of  the  intersection  in  the  local  vertical 
system.  The  ”Z"  coordinate  is  previously  known  as  the 
tidal  elevation.  These  local  coordinates  are  then 
transformed  back  into  geocentric  coordinates  to  compute  the 
constant  coefficients  of  the  collinearity  condition 
equations  in  preparation  for  return  to  the  calling  program. 


Computation  of  the  partial  derivatives  of  the 
collinearity  equations  now  takes  place.  As  discussed  in 
the  development  of  the  mathematical  model,  the  collinearity 
condition  does  not  hold  for  submerged  points.  However,  the 
partial  derivatives  of  these  equations  can  be  modified  to 
correct  for  the  refraction  of  the  image  rays  which 
invalidates  this  condition.  The  partials  are  computed  as  a 
function  of  the  image  ray  -  water  surface  intersection 
point.  In  this  manner,  the  submerged  point  can  be 
incorporated  into  the  standard  aerotriangulation  problem. 


The  partial  derivatives  for  the  submerged 
points  are  computed  in  the  local  vertical  system  and  then 
transformed  into  the  original  triangulation  system  in 
preparation  for  solution  of  the  triangulation  problem. 
These  partial  derivatives  are  described  in  detail  in 
Section  4.1  of  this  report.  At  this  point  control  returns 
to  the  calling  routine  for  continuation  of  the  normal 
program  processing. 

None 

/FRDATA/  / GNDATA/  /WATER/ 

/ COEFPT/  /INDX/  /FLAGS/ 

/DATAB3/  /AXES/  /TAPES/ 

None 
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REFERENCED: 


RTFRM4 ,  PLHXY1,  WATPNT,  MPYAB3,  ADDMT3 


4.2.2  SUBROUTINE  WATPNT 

SUBROUTINE:  WATPNT 


DATE  WRITTEN:  May  1980 

WRITTEN  BY:  L.  Peters. 


LANGUAGE/MACHINE:  FORTRAN  /  UNIVAC  1108  Exec  8 


DESCRIPTION: 


CALL  LINE: 
INPUT: 


OUTPUT: 

CALLED  BY: 


This  subroutine  computes  the  coordinates  of 
the  Image  ray/water  surface  Intersection  for  a  bathymetric 
Image  point. 


CALL  WATPNT  (  XCAM,  XGRD,  TIDE,  WAT,  ICOL,  J  ) 
Thru  the  argument  list  - 


XCAM 

XGRD 

TIDE 

ICOL 


J 

XWAT 


Camera  station  coordinates 
Ground  point  coordinates 
Height  of  water  surface 


Column  pointer  to  an  Image  In  the 
physical  record 


-  Internal  frame  identifier 


Image  ray/water  surface  intersection 
coordinates 


UNWTR 


METHOD: 


The  maximum  value  for  the  refraction  angle  is  defined  by 


max 


sin 


-1  1 


and  the  minimum  value  is  0.  A  binary  interpolation  is 
performed  to  compute  the  value  of  the  refraction  angle  for 
each  image  ray.  Knowing  the  refraction  angle,  the  incidence 


angle,  I, 

can  be 

computed  using  Snell's  Law.  The  image 

ray /water 

surface 

intersection  coordinates  can  be  computed 

from  : 

<*c  -  V  (zc  -  V  1 

*W 

*c 

«*c  -  V2  +  (Yc  -  V2)l/2 

V  _ 

<Yc  -  V  (zc  -  V  ““  1 

Yw 

Yc 

«XC  -  V2  +  «c  -  V2>1/2 

LIMITATIONS: 


None 


COMMON  BLOCKS 


/FLAGS/  /INDX/ 
OPTIONS:  None 

SUBROUTINES 

REFERENCED:  NONE 


/TAPES/ 
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4.2.3 


SUBROUTINE  RTLCGC 


SUBROUTINE: 

RTLCGC 

DATE  WRITTEN: 

June  1980 

WRITTEN  BY: 

L.  Peters 

LANGUAGE/MACHINE: 

FORTRAN  /  UNIVAC  1108  Exec  8 

DESCRIPTION: 

This  subroutine  computes  the  local  to 

geocentric  rotation 

for  each  bathymetric  point. 

CALL  LINE: 

CALL  RLTCGC 

INPUT: 

Thru  ITAPE2  - 

NREC 

-  Number  of  records  in  subrecord 

IPNTER 

Pointer  to  beginning  column  of 
subrecord 

JJ 

-  Total  number  of  columns  in  subrecord 

RECORD 

-  Ground,  image,  elimination, 

addition,  and  camera  distance 

subrecords 

OUTPUT: 

Thru  ITAPE2  - 

IWATRT 

-  Local  to  geocentric  rotation  matrix 
for  the  bathymetric  point 

CALLED  BY: 

LINK5A 

METHOD: 

See  MUSAT  Program  Documentation  for  subroutine 

RTFRM6 . 

LIMITATIONS: 

None 

COMMON  BLOCKS: 

/COUNTR/  /TAPES/ 

/FLAGS/ 

OPTIONS: 

None 

SUBROUTINES 

REFERENCED:  RTFRM6 
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4.3  DICTIONARY  OF  NEW  COMMON  BLOCKS  -  MUSAT-IIIA 

I  - 

COMMON  /INDX/  SN 

VARIABLE  TYPE  DESCRIPTION 


SN 

Real 

Index  of  refraction  of  the  water  body 

COMMON  /WATER/  IREC, 

IWATPT( 11 ,50) 

VARIABLE 

TYPE  DESRIPTION 

IREC 

Integer 

Number  of  point  records 

IWATPT(l.I) 

Integer 

Point  ID 

IWATPT(2,I) 

Integer 

Bathymetric  flag 

IWATPT(3-11, 

I) 

Integer 

Local  to  geocentric  rotation  matrix 
for  a  bathymetric  point 

6 
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MODIFICATION  OF  AS- 11AM  FRAME  SOFTWARE 
FOR  BATHYMETRIC  COMPILATION 

In  addition  to  the  capability  to  triangulate  bathymetric  points 
using  the  MUSAT  programs,  the  capability  to  compile  bathymetric  points  on  the 
AS-11AM  using  conventional  frame  imagery  was  provided.  The  section  develops 
the  mathematics  involved  in  the  computation  of  the  true  coordinates  of 
bathymetric  points. 


5.1 


COMPUTATION  OF  TRUE  COORDINATES  FOR  BATHYMETRIC  POINTS 


All  computations  for  the  correction  of  measured  model  coordinates 
are  performed  In  the  model  coordinate  system.  The  tide  elevation  for  each 
frame.  Input  In  a  geographic  coordinate  system  Is  transformed  Into  model  space 
by  applying  the  appropriate  transformation  coefficients.  We  can  now  define  the 
following  notation: 


B*l»  Bylf  BZj 

Bx  2 ,  By2»  BZ2 

-  Model  coordinates  of  the  left  and 
right  camera  stations  respectively. 

Zwj,  Zw2 

-  Model  space  height  of  water  surface 
on  left  and  right  frames 

respectively. 

n 

-  Index  of  refraction  of  water. 

(Assumed  equivalent  on  both  frames) 

Xp.  Yp.  Zp 

-  Model  coordinates  of  a 
bathymetric  point. 

measured 

Xwj ,  Ywj 

Xw2,  Yw2 

-  Model  coordinates  of  image  ray/water 
surface  intersection  on  left  and 
right  frames  respectively. 

Xj,  Yj,  Zj 

-  True  model  coordinates  of  a 
bathymetric  point. 

measured 

The  first  step  in  computing  the  true  model  coordinates  is  to 
compute  the  Image  ray/water  surface  intersection  for  each  frame.  We  have: 


Xw 


Yw 


Xw 


1 


1 


2 


BXj  +  (BZl  -  ZwL) 
Byj  +  (BZl  -  ZwL) 
Bx2  +  (BZz  -  Zw2) 


(Xp  -  Bxj) 
<BZj  -  Zp) 


(Yp  -  By  j) 
(B*j  -  Zp) 


(Xp  —  BX2) 

(Bz2  -  Zp) 


(5.1) 


(5.2) 


(5.3) 
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Next,  the  angles  of  incidence  and  refraction  of  the  image  rays  are  computed. 


Let, 

Dj  -  ((BXl  -  Xwj)2  +  (Byj  -  Yw^2)^  (5.5) 

D2  »  ((Bx2  -  Xw2)2  +  (By2  -  Yw2)2/^2  (5.6) 


The  angles  of  incidence  are,  then  : 


The  angles  on  refraction  due  to  the  water  surface  are  now  : 

[^] 

[-2^] 


(5.7) 

(5.8) 

(5.9) 

(5.10) 


The  true  coordinates  of  the  bathymetric  are  defined  by  the 
intersection  of  the  underwater  image  rays  of  a  bathymetric  point  from  each 


(5.13) 


Now, 

(Xj  ~  XWj)  -  (r j/r3)  (ZT  -  Zwp  -  r13zT  "  ri3Zwj 

(Yt  -  YWj)  -  ^r2^r3^  ^ZT  ”  ZwP  "  r23ZT  ”  r23Zwl 

-  Xw2)  -  (rj/rp  (ZT  -  Zw2)  =  rJ3ZT  -  rJ3Zw2 

«T  -  *»2>  '  (r2/rP  (ZT  -  Zw2>  '  r23ZT  -  r23Zw2 

In  matrix  form. 


1 

0 

“r13 

xT 

_XWj 

-  ijj*.. 

[ail 

1 

0 

~rl3 

YT 

X»2 

-  r13Zw2 

3i 

0 

1 

_r23 

ZT 

Yw 

1 

‘  r23Zwl 

a2 

0 

1 

“r23  _ 

Yw 

L  2 

-  r23Zw2 

La2. 

rJ3  -  tan  (XWj  -  Bxp  /  -Dj 

r23  "  tan  ei  (Ywi  -  Byj)  /  -0! 

rj3  -  tan  02  (Xw2  -  Bx2>  /  -D2 

r’3  -  tan  ©2  (Yv2  -  By^  /  -T>2 


(5.14) 

(5.15) 

(5.16) 


(5.17) 


(5.17a) 

(5.17b) 

(5.17c) 

(5.17d) 


The  normals  are, 


2 

0 

-<r13+  r13> 


0 

2 


-(r23+  r23) 


"(r13+  rl3) 

*T 

Va; 

~(r 23+  r23) 

yt 

a2+a2 

,  2 .  ,  2 ,  2 .  ,  2n 
(r13  +  r13  r23  +  r23  * 

ZT 

0 
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■ 

Let, 

0  -  -<rj3a|  +  r;3S  +  r23*2  +  ^3a2> 
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Solving  the  system,  we  have  the  true  coordinates  of  the  bathymetric  point 


6.0 


TESTING  AND  EVALUATION  OF  MODIFIED  MUSAT  PROGRAMS  AND  AS-11 


SOFTWARE 


This  section  presents  the  results  of  the  operational  testing  and 
evaluation  of  the  modified  MUSAT  aerotriangulation  program,  and  the  modified 
AS-11  compilation  software.  Due  to  time  constraints,  the  testing  of  the 
modified  MUSAT-IIIA  program  is  not  yet  complete  and  therefore  its  test  results 
are  not  included  in  this  report. 


TESTING  &  EVALUATION  -  MUSAT 

According  to  Contract  Section  4.2.2  Quality  Assurance,  the 
following  conditions  will  govern  the  Testing  &  Evaluation  of  the  modified 
software: 


"To  verify  the  validity  and  accuracy  of  the  modified  programs, 
the  contractor  shall  perform  a  photogrammetric  block 
aerotriangulation  using  reference  photography  and  associated 
numerical  data  provided  by  the  government.  A  minimum  of  10 
photographs  will  be  used  for  this  test.  Image  points  with 
known  ground  coordinates  will  be  carried  as  unknowns  in  the 
solution.  At  least  half  of  the  image  points  used  shall  be 
subsurface  points.  The  accuracy  goal  for  these  modified 
programs  will  be  to  achieve  photograph  residuals  with  a  root- 
mean-square  error  of  IS  micrometers  with  calm  sea  conditions". 


The  testing  and  evaluation  of  the  programs  was  performed  by  DMAHTC  and 
Autometric  personnel  at  DMAHTC. 

The  test  data  consists  of  a  10  photograph  strip  of  conventional 
frame  photography  of  the  Cayman  Islands.  The  strip  was  marked  and  measured  by 
DMAHTC  personnel,  and  triangulated,  without  regard  to  the  effects  of 
refraction,  using  the  production  MUSAT-IIA,  and  MUSAT-IIIA  programs.  There 


were  a  total  of  98  ground  points  measured  on  the  film.  Of  these,  12  were 
bathymetric  points.  All  of  the  bathymetric  points  were  located  on  four 
frames  of  the  10  frame  strip.  It  must  be  noted  that  the  strip,  without  the 
refraction  model  in  place,  would  not  triangulate  within  the  standard 
convergence  criteria  of  differential  rotation  and  differential  translation, 
for  either  MUSAT  program.  The  convergence  criteria  were  relaxed  by  one  order 
of  magnitude  in  order  to  allow  convergence  of  the  solution. 

There  are  several  factors  which  may  have  contributed  to  the 
failure  of  this  strip  triangulation  solution  to  converge.  The  first  has  to  do 
with  the  geodetic  data  used  to  control  the  strip.  The  control  data  provided 
was  believed  to  be  of  first  order  quality.  However,  after  analyzing  the 
behavior  of  the  control  points  during  the  triangulation  solution,  it  appears 
that  the  control  may  not  be  of  first  order  quality.  Secondly,  the  control 
points  were  not  easily  identifiable  on  the  photographs.  The  control  points 
consisted  of  monumented  horizontal  and  vertical  marks.  None  of  the  control 
was  paneled.  Therefore,  the  inherent  error  in  the  Identification  of  the 
control  points  during  the  point  mensuration  is  not  insignificant. 

The  standard  error  of  the  point  mensuration  is  not  well  known. 
Various  values  were  assigned  for  the  standard  error  of  the  measured  image 
coordinates  during  the  testing  phase.  As  suggested  previously  in  this 
report,  Autometric  has  suggested  assigning  an  image  coordinate  measurement 
standard  error  of  1.5  to  2  times  the  normal  value  to  measured  bathymetric 
points.  The  adopted  value  for  all  measured  image  coordinates  was  10  microns 
in  x  and  y. 

Notwithstanding  these  problems,  this  data  set  was  used  for  the 


testing 


The  following  guidelines  were  used  In  the  Testing  &  Evaluation  of 
the  modified  MUSAT  programs: 

1.  Run  a  standard  production  (non-ba thyme trie)  data  set  using  the 
production  version  of  MUSAT  and  the  modified  version  of  MUSAT 
exercising  all  input  and  output  options. 

2.  Run  a  bathymetric  data  set  through  the  modified  MUSAT  program 
exercising  all  Input  and  output  options. 

3.  Analyze  the  results  of  the  bathymetric  point  triangulation  as 
compared  with  the  production  program  triangulatlon  with  no 
refraction  correction. 

4.  Analyze  the  computed  bathymetric  heights  (depths)  as  compared  to 
sounding  data  obtained  using  standard  bathymetric  techniques. 
Also  compare  to  hand  calculation  of  estimated  depth  change  due  to 
water  column. 

RESULTS  OF  TESTING  -  MU SAT- I IA 


The  first  two  Items  in  the  testing  guidelines  are  designed  to 
verify  that  all  standard  program  operation  options  are  still  addressable. 
These  options  include  such  things  as:  variations  of  the  input  coordinate 
system,  pre-processed  or  unprocessed  plate  measurements,  card  output,  etc. 
All  of  the  options  were  verified. 

The  results  of  the  third  item  In  the  testing  guidelines  are 


summarized  In  Tables  6-1  through  6-10 


The  analyses  presented  here  in  table  form  are  the  standard 
statistical  reports  generated  by  a  MUSAT  program  run.  Tables  6-1  and  6-5  show 
the  analysis  of  the  adjusted  camera  station  positions  for  standard  MUSAT-II 
program  and  the  modified  MUSAT-II  program.  The  sample  size  denotes  the  number 
of  camera  stations  adjusted  during  the  program  execution.  The  next  entry.  Sum 
of  Differences,  gives  the  accumulated  sums  of  the  differences  between  the 
initial  estimates  for  the  camera  station  positions  and  the  adjusted  camera 
station  positions.  Note  the  improvement  in  the  position  analysis  for  the 
modified  program  versus  the  standard  program.  This  implies  the  the  modified 
program  more  accurately  represents  the  geometric  characteristics  of  the 
imaging  process  than  the  standard  program;  as  is  expected. 

Inspection  of  Tables  6-2  and  6-6  show  only  a  slight  improvement 
for  the  camera  attitudes  from  the  modified  program.  Again,  only  four  of  the 
frames  in  the  strip  had  water  points  imaged. 

Tables  6-3  and  6-7  show  the  ground  point  analysis  for  both 
programs.  The  ground  points  included  in  the  analysis  are  only  the  control 
points  for  the  triangulation  run.  There  were  a  total  of  7  horizontal  control 
points  and  a  total  of  35  vertical  control  points  in  the  strip.  The 
differences  again  are  the  differences  between  the  given  control  point 
coordinates  and  the  adjusted  values  of  control  point  coordinates.  Table  6-7 
shows  a  slight  improvement  in  the  Latitude  and  Longitude  components  for  the 
control  points.  As  expected,  the  greatest  improvement  occurs  in  the  Height 
components.  This  is  due  to  the  fact  that  the  refraction  effects  imparted  by 
the  water  column  above  a  point  has  its  greatest  effect  on  the  computed  height 
(depth)  of  a  point.  It  is  interesting  to  note  that  none  of  the  control  points 
showing  a  maximum  value  for  the  difference  is  a  bathymetric  point. 
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Tables  6-4  and  6-8  show  the  computed  RMS  plate  residuals  for  all 
frames  In  the  strip.  The  Tables  differentiate  between  points  which  were 
included  in  the  triangulation  solution  and  those  which  were  not.  The  most 
notable  entries  In  Table  6-4  are  the  values  for  the  y-coordinate  residuals  of 
points  on  Frames  841  to  845.  These  are  the  only  frames  in  the  strip  which 
contained  bathymetric  points.  Table  6-8  show  the  Improvements  in  the  point 
residuals  for  Frames  841  to  845  after  the  refraction  model  has  been  applied. 
The  effects  of  refraction  alone  do  not  account  for  the  large  values  for  the  y- 
coordinate  residuals  in  the  standard  program.  Table  6-4  shows  that  the  y- 
coordinate  residuals  for  all  frames  are  approximately  twice  the  values  of  the 
x-coordinate  residuals.  This  implies  some  sort  of  systematic  problem  in 
either  point  mensuration  or  in  the  application  of  image  coordinate 
refinements. 

The  first  column  of  Table  6-9  Bhows  the  values  for  the  geocentric 
coordinates  of  the  bathymetric  points  computed  by  the  standard  MUSAT-IIA 
program  before  correction  for  refraction.  The  approximate  coordinates  for  the 
bathymetric  points  can  be  computed  if  the  assumption  is  made  that  all  of  the 
corrections  to  the  computed  coordinates  occur  in  the  Z  component.  Given  the 
adjusted  camera  station  coordinates  and  ground  point  coordinates  from  the 
standard  run,  the  direction  numbers  and  angles  of  incidence  of  the  image  rays 
with  the  water  surface  are  computed.  The  coordinates  of  the  intersection  of 
the  image  ray  and  the  water  surface  are  also  computed.  The  approximate  change 
in  height  is  computed  by  applying  the  refraction  angle  to  the  direction 
numbers  and  intersecting  the  underwater  image  rays  for  point.  These  computed 
approximate  positions  are  tabulated  in  the  second  column  of  Table  6-9.  The 
actual  computed  positions  from  the  modified  MUSAT-IIA  program  are  tabulated  in 
the  third  colum  of  Table  6-9. 
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Table  6.1 


Camera  Station  Position  Analysis  Standard  MUSAT-IIA 
No  Refraction  Correction 
(Degrees,  Minutes,  Seconds) 


Latitude  Longitude  Elev.  (M)  | 


|  Sample  Size 

10 

10 

10 

|  Sum  of  Differences 

.281 

0 

0 

5.296 

104.71 

|  Mean  Value  of  Differences 

.028 

0 

0 

.530 

10.47 

|  RMS  of  Differences 

0 

0 

.140 

0 

0 

.830 

27.30 

|  Standard  Deviation  from  Mean 

0 

0 

.144 

0 

0 

.674 

26.58 

|  Maximum  Difference 

-0 

0 

.247 

0 

0 

1.543 

66.72 

|  Camera  with  Max.  Difference 

842 

843 

846 

Table  6.2 


Camera  Station  Attitude  Analysis  Standard  MUSAT-IIA 
No  Refraction  Correction 
(Degrees,  Minutes,  Seconds) 


|  X-Tilt  Y-Tilt  Heading  | 


|  Sample  Size 

10 

10 

10 

1 

|  Sum  of  Differences 

-0 

08 

46. 

908 

-1 

22 

2.182 

-0 

0 

1.033 | 

|  Mean  Value  of  Differences 

-0 

0 

52 

.691 

-0 

08 

12.218 

-0 

0 

.103 

|  RMS  of  Differences 

0 

2 

14. 

428 

0 

12 

30.847 

0 

0 

.792 

|  Standard  Deviation  from  Mean 

0 

2 

10, 

.360 

0 

09 

57.674 

0 

0 

.828 1 

|  Maximum  Difference 

-0 

3 

42. 

.153 

-0 

23 

50.092 

0 

0 

1.955 | 

|  Camera  with  Max.  Difference  JB45  846  842  | 
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Table  6.3 


Ground  Point  Analysis  Standard  MUSAT-IIA 
No  Refraction  Correction 
(Degrees,  Minutes,  Seconds) 


1  ~ 

Latitude 

Longitude 

Elev.  (M) 

"l 

|  Sample  Size 

7 

7 

35 

I 

|  Sum  of  Differences 

0  0 

.035 

-0  0 

.120 

4.84 

1 

|  Mean  Value  of  Differences 

0  0 

.005 

-0  0 

.017 

.14 

1 

|  RMS  of  Differences 

0  0 

.035 

0  0 

.050 

4.87 

1 

|  Standard  Deviation  from  Mean 

0  0 

.081 

0  0 

.050 

4.94 

1 

|  Maximum  Difference 

0  0 

.081 

-0  0 

.123 

-9.98 

1 

|  Point  with  Max.  Difference 

=======33=32=32322===23 22=3 3=2 

1103 

1215 
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n 

ti 

it 
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48 


Table  6. A 


Plate  Residuals  Standard  MUSAT-IIA 
No  Refraction  Correction  (MM) 


840 

846  | 

1  RMS  of  Trig. 

|  Points  -  x 

.004 

.005 

.004 

.003 

.006 

.008 

.008 

.009  .008 

.006 1 

1  RMS  of  Trig. 

I  Points  -  y 

.007 

.012 

.013 

.017 

.016 

.016 

.015 

.010  .009 

.006 1 

1  RMS  of  Points 

|  Not  Trig.  -  x 

.000 

.000 

.000 

.000 

.005 

.020 

.030 

.013  .000 

.000 

1  RMS  of  Non- 
1  Trig.  Pts.  -  y 

.000 

.000 

.000 

.000 

.139 

.138 

.108 

.203  .000 

.ooo| 

1  Rms  of  All 

1  Points  -  x 

.004 

.005 

.004 

.003 

.006 

.010 

.013 

.010  .008 

.006 1 

1  RMS  of  All 

1  Points  -  y 

333333333333333333 

.007 

333333 

.012 

.013 

333333 

.017 

333333 

.033 

.043 

.040 

.068  .009 

.006 1 

Frame  With 

Max.  RMS  **** 

Trig.  Pts. 


Frame  With 

Max.  RMS  Non-  **** 

Trig.  Points 


Frame  With 
Max.  RMS  All  Pts 


**** 


Table  6.5 


Camera  Station  Position  Analysis  Modified  MUSAT-IIA 
Refraction  Correction  Applied 
(Degrees,  Minutes,  Seconds) 


1 

Latitude 

Longitude 

Elev.  (M)  | 

|  Sample  Size 

10 

10 

10  | 

|  Sum  of  Differences 

-0  0 

.103 

0 

0 

4.798 

95.53  | 

|  Mean  Value  of  Differences 

-0  0 

.100 

0 

0 

.480 

9.55  | 

|  RMS  of  Differences 

0  0 

.144 

0 

0 

.805 

25.00  | 

|  Standard  Deviation  from  Mean 

0  0 

.151 

0 

0 

.682 

24.36  | 

|  Maximum  Difference 

-0  0 

.259 

0 

0 

1.623 

64.54  | 

I  Camera  with  Max.  Difference 

844 

844 

846  1 
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Table  6.6 


X-Tilt  Y-Tilt  Heading 


le  Size  10  10 

10  | 

of  Differences  -0  11  24.831  -1  17  4. 

377  -0  0  .4331 

|  Mean  Value  of  Differences  -0  1  4.483  -0  7  42.438 


|  Standard  Deviation  from  Mean  0  1  55.159 
|  Maximum  Difference  -0  3  55.337 
|  Camera  with  Max.  Difference  845 


846 


0  0  .947 | 


Table  6.7 


Ground  Point  Analysis  Modified  MUSAT-IIA 
Refraction  Correction  Applied 
(Degrees,  Minutes,  Seconds) 


Latitude  Longitude  Elev.  (M) 


|  Sample  Size 

7 

7 

35 

1 

|  Sum  of  Differences 

0 

0 

.010 

-0  0 

.111 

-  .66 

1 

|  Mean  Value  of  Differences 

0 

0 

.001 

-0  0 

.016 

-  .02 

1 

|  RMS  of  Differences 

0 

0 

.040 

0  0 

.042 

4.46 

1 

|  Standard  Deviation  from  Mean 

0 

0 

.043 

0  0 

.042 

4.52 

1 

|  Maximum  Difference 

0 

0 

.070 

-0  0 

.068 

-9.79 

1 

|  Point  with  Maximum  Difference 

1210 

1109 

1104 

7 

Table  6.8 

Plate  Residuals  Modified  MUSAT-IIA 
Refraction  Correction  Applied  (MM) 


44  845 


333333333*333333333 

1  RMS  of  Trig. 

IS3333 

*33333 

333333 

333333 

1  Points  -  x 

.004 

.005 

.004 

.002 

.004 

.003 

.004 

.007 

.006 

.004 1 

1  RMS  of  Trig. 

1  Points  -  y 

.007 

.011 

.013 

.015 

.013 

.011 

.008 

.004 

.006 

.oos| 

1  RMS  of  Points 

1  Not  Trig.  -  x 

.000 

.000 

.000 

.000 

.005 

.004 

.002 

.001 

.000 

.ooo| 

1  RMS  of  Non- 
|  Trig.  Pts.  -  y 

.000 

.000 

.000 

.000 

.141 

.101 

.002 

.002 

.000 

,ooo| 

1  Rns  of  All 
|  Points  -  x 

.004 

.005 

.004 

.002 

.004 

.003 

.004 

.007 

.006 

.004 1 

1  RMS  of  All 
|  Points  -  y 

.007 

.011 

.013 

.015 

.032 

.031 

.008 

.004 

.006 

.oos| 

Praise  With 

Max.  RMS 

Trig.  Pts. 

**** 

Praise  With 

Max.  RMS  Non- 
Trig.  Points 

HA* 

Table  6.9 


Comparison  of  Computed  Bathymetric 
Geocentric  Point  Coordinates  (H) 


3=SSSSS33S=3 

I  Point 

3SS333333 

333333S3S3333335 

MUSAT-IIA 

Predicted 

MUSAT-IIA 

j  Number 

No-Ref rac. 

Pos.  Approx. 

Refraction 

X 

907446.7 

907446.7 

907450.4 

1047* 

Y 

-5951412.7 

-5951412.5 

-5951413.1 

Z 

2099645.1 

2099643.9 

2099644.8 

X 

904464.1 

904464.1 

904464.7 

1056* 

Y 

-5952644.9 

-5950278.3 

-5952644.5 

Z 

2097470.7 

2097470.2 

2097466.9 

X 

911407.6 

911405.7 

911413.8 

1042 

Y 

-5950277.0 

-5950277.4 

-5950279.9 

Z 

2101164.5 

2101162.6 

2101164.7 

X 

911393.6 

911393.9 

911400.0 

1043 

Y 

-5949788.3 

-5949788.4 

-5949795.4 

Z 

2102554.4 

2102558.8 

2102555.2 

X 

910095.7 

910096.4 

910101.4 

1044 

Y 

-5949944.2 

-5949944.9 

-5949940.9 

Z 

2102655.6 

2102663.6 

2102656.6 

X 

909448.2 

909448.2 

909453.9 

1046 

Y 

-5951347.2 

-5951347.2 

-5951346.9 

Z 

2098993.3 

2098992.4 

2098992.1 

X 

907741.1 

907740.2 

907745.3 

1049 

Y 

-5950406.7 

-5950406.7 

-5950409.4 

Z 

2102371.8 

2102382.8 

2102374.6 

X 

907876.0 

907878.2 

907880.4 

1050 

Y 

-5950185.4 

-5950187.2 

-5950190.0 

Z 

2102946.5 

2102963.9 

2102950.1 

X 

905938.3 

905938.9 

905940.3 

1051 

Y 

-5950307.8 

-5950308.2 

-5950302.6 

Z 

2103416.2 

2103422.9 

2103418.7 

X 

904873.0 

904871.9 

904874.1 

1052 

Y 

-5950709.9 

-5950709.8 

-5950706.3 

Z 

2102743.1 

2102722.8 

2102745.1 

1 
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Table  6.9  (cont.) 


ri 

4  1 


I 


X  904824.7  904819.1  904826.4 
Y  -5950463.2  -5950459.6  -5950462.7 
Z  2103484.2  2103490.2  2103488.2 


X  903378.8  903379.0  903378.3 
Y  -5950270.4  -5950270.7  -5950272.0 
Z  2104631.8  2104636.4  2104636.3 


Point  not  triangulated. 

Point  is  a  vertical  control  point. 


p 


i 


C 
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TESTING  &  EVALUATION  -  AS-11  SOFTWARE 


According  to  Contract  Paragraph  F.2  (Modified): 

"Acceptance  or  rejection  of  the  software  program  shall  be 
determined  using  the  following  criterion.  The  positions  of 
points  obtained  using  the  software  program  developed  under  this 
contract  modification  shall  not  differ  significantly  from  the 
positions  obtained  for  the  same  points  using  the  modified  MUSAT 
software  program  developed  under  this  contract.  The  contours 
obtained  by  the  AS-11AM  using  the  delivered  software  program 
shall  compare  favorably  with  those  of  the  same  area  obtained  by 
standard  bathymetric  means. 

The  test  data  consisted  of  a  three  photograph  strip  of 
conventional  frame  photography  of  St.  Thomas,  USVI.  The  strip  was 
triangulated  using  the  modified  MUSAT-ILA  program.  The  resulting 
triangulation  data  was  used  as  input  to  the  modified  FCA  Frame  Compilation 
software  on  the  AS-11.  The  following  guidelines  were  used  in  the  Testing  & 
Evaluation  of  the  modified  AS-11  compilation  software: 

1.  Obtain  MUSAT-IIA  or  MOSAT-IIIA  triangulation  data  and  the 
accompanying  film  for  the  stereomodels. 

2.  Perform  a  stereomodel  setup  using  the  results  of  the  MUSAT 
triangulation,  and  verify  all  standard  FCA  options. 

3.  Collect  geomorphlc  (point)  and  profile  data. 

4.  Analyse  the  collected  geomorphlc  data  as  compared  to  the 
triangulation  results  of  the  MUSAT  run. 

5.  Analyse  the  collected  profile  data  as  compared  to  sounding  data 


collected  using  standard  techniques  (NOS  sounding  sheets). 

RESULTS  OF  TESTING  -  AS-11  SOFTWARE 

The  St.  Thomas  triangulation  data  from  the  MUSAT  run  was  used  as 
Input  to  the  modified  AS-11  program.  A  model  setup  was  performed.  After  the 
Interior  orientation  was  performed,  the  stereomodel  set  up  satisfactorily, 
that  Is  with  a  minimum  of  y-parallax.  Photo  coordinates  on  the  AS-11  were 
recoverable  to  within  10-15  microns  of  their  final  values  from  the  MUSAT  run. 

The  comparison  of  the  geomorphic  data  collection  and  the  MUSAT 
triangulation  shows  differences  in  the  position  vectors  of  0.5  -  1.6  meters; 
with  the  greatest  difference  being  in  the  Z  component.  Some  of  the  factors 
affecting  the  ability  to  perceive  stereo  underwater  are:  "shimmer"  of  the 
surface  due  to  unfavorable  sun  angle,  reflectance  characteristics  of  the 
bottom,  and  the  inability  of  the  cartographer  to  recognize  when  x-parallax  has 
been  completely  removed,  i.e  when  the  "dot"  is  on  the  bottom. 
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